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A general formulation of scalar hysteresis is proposed. This formulation is based on two steps. First, a generating 
function g(x) is associated with an individual system, and a hysteresis evolution operator is defined by an appropriate 
envelope construction applied to g{x), inspired by the overdamped dynamics of systems evolving in multistable free 
energy landscapes. Second, the average hysteresis response of an ensemble of such systems is expressed as a functional 
integral over the space G of all admissible generating functions, under the assumption that an appropriate measure has 
been introduced in G. The consequences of the formulation are analyzed in detail in the case where the measure is 
generated by a continuous, Markovian stochastic process. The calculation of the hysteresis properties of the ensemble is 
reduced to the solution of the level-crossing problem for the stochastic process. In particular, it is shown that, when the 
process is translationally invariant (homogeneous), the ensuing hysteresis properties can be exactly described by the 
Preisach model of hysteresis, and the associated Preisach distribution is expressed in closed analytic form in terms of 
the drift and diffusion parameters of the Markovian process. Possible applications of the formulation are suggested, 
concerning the interpretation of magnetic hysteresis due to domain wall motion in quenched-in disorder, and the 
interpretation of critical state models of superconducting hysteresis. 



PACS: 75.60.Ej 02.50.Ga 05.40.+j 

I. INTRODUCTION 

The study of hysteresis has been a challenge to physicists and mathematicians for a long time, hi physics, 
hysteresis brings all the conceptual difficulties of out-of-equilibrium thermodynamics [1-5], first of all the 
fact that we do not know the general principles controlling the balance between stored and dissipated energy 
in hysteretic transformations [6,7]. In mathematics, on the other hand, the central issue is the formulation of 
sufficiently general mathematical descriptions grasping the essence of hysteresis beyond the limited interest 
of ad hoc models [8-11]. 
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In this article, we introduce and discuss a formulation of hysteresis of some generality, inspired by the 
following situation, often encountered in physical systems. We know that in physics hysteresis is the 
consequence of the existence of multiple metastable states in the system free energy F(X) (the temperature 
dependence is tacitly understood), and of the fact that the system may be trapped in individual metastable 
states for long times. Let us consider the simple case where the state variable Z is a scalar quantity and the 
relevant free energy in the presence of the external field H is G{X;H) = F{X) - HX. The metastable states 
available to the system are represented by G minima with respect to X, for which dO/dX = 0, ^G/dX^ > 0. 
When H is changed over time, the number and the properties of these minima are modified by the variation 
of the term -HX. The consequence is that previously stable states are made unstable by the field action and 
the system moves to other metastable states through a sequence of Barkhausen jumps. Because the condition 
dG/dX = is equivalent to H= dF/dX, one can analyze the problem by using the field representation shown 
in Fig. 1 . The response of the system, expressed in terms of H(X), is obtained by traversing the upper and 
lower envelopes to dF/dX shown in the figure, the former and the latter applying to increasing and decreasing 
H, respectively. From the physical viewpoint, this construction amounts to assuming that the system, once 
made unstable by the action of the external field, jumps to the nearest available energy minimum, which 
means that one excludes the presence of inertial effects, that could aid the system to reach more distant 
minima. 

The method discussed in Section n translates this picture into a well-defined mathematical formulation, 
based on the following two steps. 

(i) Given the time-dependent input ht and the continuous function g(x), analogous to the free energy 
gradient dF/dX of Fig. 1, one associates with them a certain evolution operator T[g](ht), which 
expresses in mathematical terms the envelope construction of Fig. 1 (Section n.l). The function g 
will be called the generating function of T[g](h,). The evolution operator acts on a given initial state Sq 
associated with the initial input ho and transforms it into the final state s, = T[g]{h,) Sq. 

(ii) Let G represent the functional space of all admissible generating functions. Then one constructs a 
general hysteresis operator as a parallel connection of the collection of operators Tig](h,) obtained by 
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varying g over G, with appropriate weights described by some measure /J. on G. Hence, one arrives at 
the following formulation, in which the overall state S, describing the collection is expressed in terms 
of the functional integral 



where sq may depend itself on g (Section n.2). 

The generality of the formulation comes from the general nature of the space G as well as from general 
ways of assigning a measure on this space. We will show that several known mathematical descriptions of 
hysteresis, like the Preisach model [9,12], are particular cases of Eq.(l), and we will discuss some new 
connections that emerge from the broader perspective offered by the functional integral formulation. A case 
of particular interest to physics is when Eq. (1) is interpreted as the average hysteresis response of a 
statistical ensemble of independent systems, each evolving in a different free energy landscape. The space G 
acts then as a probabihty space and the measure n describes the probability that an individual system of the 
ensemble is characterized by a particular generating function g e G.ln this case, there are situations that can 
be analytically investigated to a considerable degree of detail, first of all the one where the generating 
functions g(x) are interpreted as sample functions of a continuous Markovian stochastic process (Section III). 
In particular, we will show that homogeneous processes give rise to Preisach-type hysteresis, and we will 
derive explicit analytical expressions for the Preisach distribution as a function of the parameters governing 
the statistics of the Markovian process. 

The results obtained in this article can be important in apphcations to physics, where randomness due to 
structural disorder often plays a key role in the appearance of hysteresis effects. The equivalence between 
Markovian disorder and Preisach-type hysteresis implies that the average system response under small fields 
is parabolic, a result well known in magnetism under the name of Rayleigh law [12]. In superconducting 
hysteresis [13], the same equivalence might be of help in the interpretation of critical state models [14,15] in 
terms of the statistics of the pinning sources acting on Abrikosov vortices, given the equivalence between 
this class of models and the Preisach model [9,16,17]. Conversely, a hmitation of our formulation is the fact 
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that it is based on independent single-degree-of-freedom subsystems, and is thus expected to yield an 
incomplete description of hysteresis effects arising in systems with more complex internal structures 
[7,18,19]. 



II. MATHEMATICAL DESCRIPTION OF THE MODEL 

In this section, we will discuss in more detail the various ingredients defining the structure of the model: 
input histories, generating functions, admissible states, input-output relationships, stability properties, and 
finally the functional integration over G. 

II.l Hysteresis in individual systems 

Let us consider an individual system characterized by a particular function g{x). The system is acted on 
by the scalar time-dependent input ht and generates the scalar output x,, in a way dependent on the function 
8(x). 

A. Input histories. We shall consider input histories h(t), t>0, such that, at any time, < h(t) < hy, where 
and hu are fixed given fields, delimiting the input range of interest. They will be termed lower and upper 
saturation field, respectively. The function h(t) will be assumed to be piece -wise monotone. 

B. Generating functions. Let us consider a given output interval [xz,,X[/]. We say that the function g{x) is an 
admissible generating function associated with the interval fex;/] if it satisfies the following properties (see 
Fig. 2): 

(i) g is continuous in [j:z,,X£/] 

(ii) g(xL) = Hl, g(xu) = hu (2) 

(iii) hL < g(x) < hu for any x in the range Xl<x<Xu 

C. States. Any ordered input-output pair s - (h,x) is an admissible state for the system. We will be mainly 
interested in equilibrium states, defined as the states of the form s = (h = g(x),x), with Xl < x < Xu. In other 
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words, an equilibrium state is represented by a point on the generating function. All other states will be 
generically termed jump states. Given the generating function g(x), an equilibrium state is fully described by 
its output X. In this sense, we will often identify an equihbrium state (h = g(x),x) simply by its x value. When 
the input h is given, the possible equilibrium states under that input are obtained by solving the equation g{x) 
= h. In general, more than one solution will exist. Only the states Sl - {hiyXi) and Sv - (hu,Xu) are unique by 
definition. They will be termed lower and upper saturation state, respectively. We will assume that the state 
of the system before any action is made on it is always an equilibrium state. 

D. Auxiliary functions. Given the generating function g{x) and the equilibrium state Xq, let us introduce the 
function h[g](x;xo), defined as (see Fig. 2): 



min g if Xj^<x<Xq 

(3) 

max g if Xq<x<Xjj 



where the symbols "min" and "max" indicate the minimum and the maximum of g{x) in the specified 
interval. Function h[g] has the character of a non-decreasing envelope to g(x), more precisely, of an upper 
envelope for x>Xo and a lower envelope for x < Xq. The inverse of h[g]{x;xo) will be denoted by X[g](/i;xo). The 
mathematical aspects of the connection between h[g] and X[g] are discussed in Appendix A. 

E. System evolution. Let us introduce the following evolution operator T[^](/i), defined in terms of the 
function X[g](h;xoy. given the equihbrium state Sq = (ho,xo), with ho = g(xo), and the input value h, the state s 
obtained by applying the input h to So is given by the expression 

^ = ■^o = (h, -^[g] (h; Xq )) (4) 

The evolution of the system is constructed by applying T[g](h) many times in sequentially, once for each 
input reversal, as shown in Fig. 3. More precisely, let us suppose that at the initial time f = the system is in 



5 



the equilibrium state Sq = {ho,Xo) and let us apply the piece-wise monotone input history h{t). Let us denote by 
hi, h2, ...,h„ the sequence of input values at which the input is reversed in the time interval [0,t], and finally 
let h, be the current input at the time t. Then, the state of the system at the same time is given by 



(5) 



hi particular, the output x, can be expressed in the form 



(6) 



where x„ - the output value at the last reversal point - depends in general on all past reversal inputs. Note that 
the evolution is rate-independent, because the state s, depends only on the current value of the input and on 
the sequence of past reversal inputs, regardless of the input rate of change. 

A relevant aspect of the evolution described by Eq. (5) is that it exhibits return-point memory (also called 
wiping-out property) [19-21]. By this we mean that, given the initial equilibrium state So = (haXo) and the 
input extrema sequence ho, hi, hi, hi, with hi > ho, ho< h2< hi, then T[g](/zi) T[g](h2) T[g](hi) so = T[g](hi) so 
(identical conclusions apply to the case hi < ho, ho> h2> hi). In other words, when the input returns back to 
the first reversal value hi, the system returns back to the exact same state it occupied when the input first 
reached the value hi, and the effect of the intermediate input extrema is wiped out. To prove the existence of 
return-point memory, we begin by remarking that return-point memory is a property of any system whose 
time evolution satisfies the following properties [19,12]: 

(i) the evolution is rate-independent; 

(ii) there exists a (partial) ordering relation among the states of the system; 

(iii) ordering is preserved during the evolution of the system under the action of ordered input histories. 



Property (i) is the direct consequence of the definition of the evolution operator T[g](/?). As property (ii) is 
concerned, there exists a natural ordering relation deriving from the fact that an equihbrium state is identified 
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by its output value. In fact, given the equilibrium states Si = (hi^i) and S2 = {h2,X2), we can simply state that 
si < S2 if xi < X2 in the usual sense. Finally, property (iii) is the consequence of the theorems of Appendix B, 
which show that the ordering just defined is preserved by Eq. (5) under the application of ordered input 
histories. Therefore, return-point memory is indeed a property of Eq. (5). 

As discussed in [9, p. 13], return-point memory has the consequence that the final state s, defined by Eq. 
(5) is controlled (assuming, for simplicity, that the system is initially in the lower saturation state) by the 
alternating sequence of dominant extrema Hmi, h^i, hm, h^i, ... contained in the full reversal sequence hy, /i2, 
... , /i„. By this we mean that hmi is the global input maximum in the time interval [Q,t\, hm\ is the global 
input minimum in the time interval [tM\,i\, where tui is the time at which hui is reached, and so on. 

F. Stability properties. They can be conveniently described by introducing the concepts of strong and weak 
stability. Given the equilibrium states Sa - (/za,xa) and Sb = {hB,x^, with xa < Xb, we will say that Xa and Xb 
belong to the same strongly stable (strongly unstable) interval if Ha < hs (Ha > hs) and T[g](h) Sa = T[g](h) Sb 
for any input h. This definition introduces an equivalence relation among equilibrium states, which generates 
a partition of the interval [xl,Xu\ into disjoint strongly stable and strongly unstable sub-intervals, possibly 
reducing to isolated points. Conversely, given the equilibrium states Sa = (^a.Xa) and Sb = (hB,XB), with Xa < Xb 
and Ha < hs, we will say that the segment [a:a,Xb] is weakly stable if gfe) = Ha, gfe) = ha, and Ha < g{x) < 
JtB for any x in the range Xa<x < Xb. Notice that the interval [x^Xu] is always weakly stable by definition. 
When the equilibrium states Sa = (HajXa) and Sb = {hB,x^ are such that Xa < Xb and Ha > Hb, we will call the 
segment [xa,jcb] weakly unstable if gfe) = Ha, gfe) = Iib, and Ha > g(x) > hs for any x in the range Xa<x< Xb. 
Each weakly stable or weakly unstable interval will contain in general several strongly stable and strongly 
unstable sub-intervals. The various possibilities are shown in Fig. 4. One can verify from Fig. 3 that, given 
any two subsequent reversal points {h^, X),) and (hk+i, x^+i) associated with a certain input history, the output 
interval [x^,x*^+i] is always weakly stable. In a sense, the evolution operator T[g](/?) provides a mechaiusm to 
select the weakly stable portions of the given generating function. This feature will play an important role in 
the general formulation of Section 11.2 and in the particular cases discussed in Section III. 

Stability considerations are important, because the evolution of the system under varying h is reversible 
inside each strongly stable sub-interval, so that its hysteresis properties are essentially governed by the 
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sequence of jumps occurring from one stable sub-interval to another. A system initially occupying a state 
inside a weakly unstable interval will never be able to come back to this interval if it ever leaves it. 
Therefore, only the strongly stable sub-intervals that do not belong to any weakly unstable portion of fex;/] 
control the permanent hysteresis properties of the system. Two generating functions possessing the same set 
of weakly unstable sub-intervals and differing only in their values inside these intervals will give rise to 
identical hysteresis properties. Considerations of this kind permit one to recognize certain quahtative aspects 
of hysteresis independent of the details of g. For example, Preisach-type hysteresis, briefly discussed in the 
next sub-section, arises from generating functions containing two strongly stable intervals separated by a 
weakly unstable part [22]. 

II.2 Hysteresis in system ensembles 

Let us now consider an ensemble of systems of the type discussed in the previous sub-section. Each 
system is identified by a particular generating function g(x) whose domain [xz^x^] will be in general different 
from system to system. We wish to investigate the global hysteresis properties that we obtain when we 
subject the individual systems of the ensemble to the common input ht and, roughly speaking, we sum up 
their responses. The formalism whereby we will carry out this sum in precise mathematical terms is the 
following. Let us suppose that the response of each individual system is described by some quantity q[g], 
dependent on the generating function g. The ensemble value of that same quantity, say Q (we will use capital 
letters to denote ensemble properties), will be expressed as a functional integral of the form 



Equation (7) is to be interpreted in the following way. The symbol G denotes the functional space of all 
generating functions satisfying Eq. (2) for some Xl and Xu- hi general, Xl and Xu will be different for each g s 
G. One introduces a convenient set (a so-called <7-algebra) of subsets A c G and a positive measure defined 
over that algebra, fi{A) > 0. Then one assumes that there exist elements of the algebra giving rise to values of 
q[g\ inside any arbitrarily small neighborhood of a given value q[g\ = x, and uses the measure of these subsets 
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to calculate the Lebesgue integral of q[g] over G, represented by Eq. (7). To make this loose description 
mathematically rigorous, one should resort to the language and the methods of measure theory [23,24]. 
However, it is not the purpose of this paper to go deeper into these mathematical aspects. In the following 
analysis, it will be sufficient to assume that Eq. (7) does have a precise meaning as a functional integral, and 
that one knows how to assign the measure /i in specific cases. In Section III, we will discuss a particular case 
where one explicitly constructs the measure /J. and expresses the result of the functional integration in a 
closed analytic form. 

As a first step, let us apply Eq. (7) to the definition of ensemble equihbrium states. The main difference 
with respect to Section 11.1 is that we can no longer identify an equilibrium state by its output value. In fact, 
given the individual output x, the corresponding input g(x) may not exist for certain g functions (if x is 
outside the function domain [xuXu]), or may be different from function to function, which is not compatible 
with the assumption that the entire ensemble is driven by a common input history. In fact, in order to 
construct a meaningful equilibrium state we must: (i) specify the input value ho, (ii) determine, for each 
generating function g s G, the set of solutions of the equation g(x) = ho; (iii) for each g, select one of these 
solutions, say ^[g]{ho), according to some rule, and build the state S[g]{ho) = (ho, ^[g]{ho)); (iv) construct the 
ensemble state So as 



Equation (8) shows that a great number of possible equilibrium states are associated with a given input ho, as 
a consequence of the various possible choices for ^[g]{ho). Only the lower and upper saturation states are 
unique, because the equations g(x) = hi and g(x) = hu admit just one solution, xi and Xu, for each g. 
The ensemble evolution is obtained by applying Eq. (7) to Eqs. (4)-(6), that is. 




(8) 
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Because return-point memory is a property of each individual system, it will also be a property of the 
ensemble evolution. 

The formulation summarized by Eqs. (9)-(ll) is rather general and powerful, but it is also quite abstract. 
It is not obvious how one could possibly determine the measure associated with particular cases and carry 
out the functional integrals. In this cormection, a situation of interest is when one is deahng with a statistical 
ensemble of independent systems, and one wishes to calculate statistical averages over the ensemble. In that 
case, Eq. (7) translates into mathematical terms the physical idea that Q represents the sum of all the 
individual contributions q^g], each weighed by its probability d^(g) to occur. Accordingly, G must be 
endowed with the structure of a probabihty space: the elements A of the (7-algebra represent the admissible 
events that may occur in experiments, the measure ^ satisfies the postulates of probability, and ^(A) 
represents the probabihty of the event A. 

Probabihty considerations permit one to express Eqs. (9)-(ll) in the following useful form. Let us 
consider for simplicity the case where the ensemble is initially in the lower saturation state. This ehminates 
from all equations the complicated dependence on the iiutial state |[g](^o) of the individual systems. In 
particular, Eq. (11) can be written as 



At the end of the previous section, we mentioned the fact that the input reversal sequence hi, hi, ... , h„ 
selects a sequence of weakly stable portions of the generating function. Let us denote by p(x„h,;h„,...,hi) dx, 
the probability of having a function g e G such that g{u) = h, for some u in the interval [x„x + dx,] and such 
that there exists a sequence of x values, xi, X2, ... , x„, for which g(xi) = hi, g(x2) = hi ... , g(x„) - h„ and 
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[jCi^i], [xi,X2\,..., [x„.i^„], [x„^,] are all weakly stable sub-intervals. Then one can formally write Eq. (12) in 
the equivalent form 



= j-^/ p(x,,h/\,...,hi) dx, 



(13) 



The interest of Eq. (13) lies in the fact, discussed in Section HI, that the probability density p can be 
explicitly calculated in the case where the measure is generated by a continuous Markovian stochastic 
process. 

We conclude this section by showing, as an example, when Eqs. (9)-(ll) can contain and reproduce other 
known hysteresis models. We will discuss the Preisach model [9,12]. To this end, let us consider the case 
where the integral of Eq. (7) is restricted to the subspace GpczG containing the generating functions of the 
type shown in Fig. 5. The domain [xLyXu\ is equal to [-1,1] for all functions. Each function is made up of two 
strongly stable, vertical branches [25], separated by a central, weakly unstable interval. The left branch 
increases from h = h^io h = a ^.i x = Xl = and the right one increases from h = b io h = hu x = Xu = \. 
One must assume a > ^ if the central part is to be weakly unstable. Then, let us decompose the space Gp into 
the equivalence classes Aab containing all the generating functions characterized by the same a and b, and let 
us express Eq. (12) as an integral over those equivalence classes, that is 



As discussed at the end of Section 11.1, all generating functions characterized by the same set of weakly 
unstable intervals and differing only in the values they take inside these intervals give rise to identical 
hysteresis properties. This means that the function X[g} appearing in Eq. (14) takes the same values for any g 
G Aab, so it can be taken out of the integral. We obtain 




(14) 
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a>b 



(15) 



where YaAhA expresses in simplified operator form the dependence of on a, b, and input history, whereas 
H{a,b) represents the measure of the class Aj^,. It is easy to check through Figs. 3 and 5 that Yablh,] = ±1. In 
other words, Yabiht] is a rectangular-loop operator and the hysteresis model is a weighted superposition of 
these operators, which means that it is precisely the Preisach model. 

In the next section, we will show that the Preisach model can also emerge in a completely different 
context, when the generating functions g{x) are interpreted as sample functions of a Markovian stochastic 
process. 

III. MEASURES GENERATED BY STOCHASTIC PROCESSES 

The main difficulty of the formulation discussed in the previous section hes in its abstract nature. One 
needs some tools to generate and manipulate the measure before one can apply the approach to specific 
situations of interest. In this section, we discuss the case where this issue is addressed by interpreting the 
generating functions g{x) as sample functions of some stochastic process. We will show that, quite 
remarkably, the calculation of Eq.(13) can then be reduced to the solution of the level-crossing problem (also 
called exit problem or first-passage-time problem) [26,27] for the stochastic process considered. This will 
create a direct bridge between two such distant fields as the theory of hysteresis and the theory of stochastic 
processes, and will permit us to exploit the machinery of level-crossing analysis to derive analytical results 
on hysteresis. In particular, we will show that homogeneous continuous Markovian processes give rise to 
Preisach-type hysteresis and we will derive exphcit analytical expressions for the associated Preisach 
distribution n(a,b). 

III.l. Markovian processes with continuous sample functions 

Let us consider the stochastic process gj^. To avoid confusion, we point out that the independent variable x 
has nothing to do with the real time t: it will play the role of a fictitious time to be eventually identified with 
the system output. We assume that the process is Markovian, that is, its evolution under given initial 
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conditions, say atx = Xq, depends on these conditions only and not on the behavior of the process for x < Xq. 
In addition, we assume that the process is a diffusion one, which means that (almost) all its sample functions 
g(x) are continuous functions of x. We will use the letter h, with appropriate subscripts, to denote values 
taken by these sample functions. 

In Section 11.1.F, we discussed the fact that, given the generating function g(x), any arbitrary input 
reversal sequence hi, hi, ... , h„ selects a sequence of weakly stable portions of that function. We will show 
now that, when g(x) is interpreted as a sample function of g^, weakly stable intervals are naturally and 
intimately related to the solution of the level-crossing problem for g^. To this purpose, let us consider the 
interval [b, a] of the h axis and let us select in it the point ho, with b <ho<a. Let us imagine that we generate 
a sample function g(x) of the process starting from (ho,xo), and that we follow it until it reaches one of the 
two boundaries, h = b or h = a, for the first time (Fig. 6a). The value of x at which g{x) reaches the boundary 
is a random variable. The problem of determining the statistical properties of this random variable is known 
in the hterature as level-crossing problem or exit problem or first-passage-time problem for the stochastic 
process [26,27]. Let us restrict the level-crossing analysis to the sample functions that reach the upper 
boundary h = a first, and let us take the limit ho b, as shown in Fig.6b. If we interpret the function shown 
in Fig. 6b as a portion of some generating function extending outside the interval [j:i„xj, we immediately 
recognize that the interval {xi„Xa\ is weakly stable (see Section II.l.F), because g{xb) - b, g(xj = a, and b < 
g{x) < a for any x in the range Xb<x < Xa. Therefore, Xa and Xb are admissible reversal outputs that may be 
encountered under input histories with input reversals ath = a and h = b, and the solution of the particular 
level-crossing problem shown in Fig. 6b is accordingly expected to give direct information about the 
probability distribution of those reversal outputs. 

To analyze in detail the consequences of this idea, let us assume that the particular level-crossing problem 
of Fig. 6b has been solved, so that we know the conditional probability density T{a,Xa I b,Xb) of having a 
level-crossing event atx = Xa (that is, of having g(Xa) = a) conditioned by the fact that g{Xb) = b. The function 
T is defined for a>b and is non-anticipating, that is, T(a,Xa I b,Xb) = for Xa < Xb. It obeys the normalization 
condition 
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^T(a,x^ \b,Xi,) dx^ =1 



(16) 



The quantities Xa and xi, are in general random variables. Let us denote by PaiXa) and Pbixb) their probabihty 
distributions. These distributions are not independent, because they must satisfy the equation 



Notice that, because of the Markovian character of the process, Eq. (17) is fully independent of the behavior 
of the process outside the interval [jc^Xa]. Let us define the space G of Section n.2 as the space containing all 
those sample functions of the given Markovian process which satisfy the requirements of Eq. (2) for some xl 
and xu, that is, gixj) = hi, g{xu) = hu, and hi < g{x) < hu for any x in the range Xl<x<Xu-^ general, Xl and 
Xu will be random variables, taking different values for each g s G. Then, let us study the evolution of the 
ensemble described by Eq. (12), assuming that the ensemble is initially in the lower saturation state (i.e., h(0) 
= hj). Let us denote by h\, hi, ... , hn the alternating sequence of dominant input extrema controlling the 
evolution of the ensemble (this sequence was indicated as Hmu hmu Hmi, hmi, ... in Section 11.1.E). Let us 
analyze in some detail what happens when the input ht increases from Hl up to h\ along the first hysteresis 
branch, and then decreases from h\ to hi along the second one. We denote by x, the output value associated 
with h, for a given generating function (see Fig. 7a). The interval [xl,xJ is weakly stable for each g e G, so 
we can apply Eq. (17), with b = hu Xb = XL,a = h„ Xa = X[g]{h,;xL) = x„ pbiXb) = Pl(xl), Pa(Xa) = pixM- 



The probabihty distribution Pl(xl) of the lower saturation output Xl can be chosen at will, it is part of the 
characterization of the initial state of the ensemble. After that, Eq. (18) permits one to calculate the unknown 
distribution p{x„h,) on the basis of the known functions Pl{xl) and T{h,,x, I h^Xi). The distribution p{x„ht) is 




(17) 



p(x,,h,)= ^T(h,,x,\h^,Xi^)p^(xJ dx^ , h^<h,<hi 



(18) 
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exactly the function needed in Eq. (13) to calculate the average response of the system, according to the 
expression 



Equation (18) can be also used to calculate the probability distribution pi{xuh{) of the reversal output Xi at hi. 
We find 



Similar considerations apply to the second, decreasing input branch, where h, decreases from hi to hi and x, 
accordingly decreases from xi to xi (see Fig. 7b). The weakly stable interval to consider is now [xtM^- By 
applying Eq. (17) to this interval, one obtains 



The main difference with respect to Eq. (18) is that the unknown distribution p(x„h,;hi) is now inside the 
integral on the right-hand side of Eq. (21), so Eq. (21) is actually an integral equation for p{x„ht,hi). It is this 
difference in the structure of Eqs. (18) and (21) that is responsible for the onset of hysteresis in the average 
output. The comparison of Fig. 7a with Fig. 7b gives a pictorial illustration of this difference. Although the 
the probability distributions of Xl and Xi are the same, the level-crossing problems to solve under increasing 
or decreasing input are different, and therefore give rise to different probability distributions and different 
average outputs. 

The procedure that we have described can be continued to calculate the distribution piixiMiM) of the 
second reversal output, given by the solution of the integral equation 




(19) 




(20) 




(21) 
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(22) 



and then the distributions p{x„h,;h2,hi), Piix^ihiMiMi) and so on up to the distribution p„(x„;h„,...,hi) of the 
last reversal output. At this point, the probability density p{x„h,; h„, ... , hi) of the current output at time t is 
given - depending on whether the current input is increasing or decreasing - by one of the two equations, 

p(x,,h,;h„,...,h,)= ^T(h,,x,\h„,xjp„(x„;h„,...,h,) dx„ , h„<h,<h„_, 

(23) 

Pn(Xn'A'-A)= \T(h^,xJ\,x^)p(x^,h/,h^,...,h,) dx, , \_^<h^<\ 

and the corresponding average output is 

X^{h^-\,...,1\)= |x, p(x,,h,;h„,...,hi) dx, (24) 

By the analysis just concluded, we have reduced the original functional integral over G (Eq. (12)) to a chain 
of integrals and integral equations (Eqs. (20), (22), (23)), dependent on the saturation distribution Pl(xi) 
(arbitrarily chosen) and the transition density T{a,Xa I b,Xb). The central problem is then the calculation of 
T(a,Xa I b,Xb) for a given process. 

III.2. Homogeneous processes 

Particularly simple and interesting results are obtained when the statistical properties of the process 
considered are translationally invariant with respect to x, that is, when the process is homogeneous with 
respect to x. In fact, in this case it is not necessary to determine the complete function T(a,Xa I b,Xb) in order to 
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predict the hysteresis properties of the ensemble. To clarify this point, let us come back to the first of Eqs. 
(23). Because of the assumed homogeneity of the process, T(a,Xa I b,Xb) = T{a,Xa-Xb I b,0). Therefore 

X 

p{x,h-\,...,}\)= ^T{h,x-x^\h^,Q)p„{x^-\,...,}\) dx^ (25) 

where we have dropped for simplicity the t subscript in x and h. According to Eq. (24), the average system 
response is obtained by multiplying both members of Eq. (25) by x and by integrating over x. By expressing 
xs&x = {x-x,^+x„ and by rearranging the appropriate integrals on the right-hand side, we obtain 

X{h-\,...,\) = X^+X^{h\hJ ,h>h„ (26) 

where 



Xj.{a\b)= T{a,u\b,0)du , a>b (27) 



and 



= Pn{^n'\^-A)dx^ (28) 



When the second of Eqs. (23) is the relevant equation, by perfectly similar considerations one obtains 



X{h-\,...X) = ^n-^AK\h) , h<h„ (29) 



We see that the hysteresis properties of the system are fully controlled by the first moment of T only, given 
byEq. (27). 
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Equation (26) shows that the shape of a generic ascending hysteresis branch starting from the reversal 
field h„ is the same regardless of the past input history. The influence of past history is summarized in the 
value of X„ (Eq. (28)), and the branches generated by different histories differ by a mere shift along the X 
axis. The same is true for descending branches (Eq. (29)). The importance of this result lies in the fact that it 
implies the validity of the so-called congruency property [9]. It is known that return-point memory (built in 
the description from the beginning, see Section n.l E) and congruency represent the necessary and sufficient 
conditions for the description of a given hysteretic system by the Preisach model [27,9,12]. Therefore, we 
conclude that the hysteresis generated by a homogeneous, diffusion Markovian process is of Preisach type. 
The process is fully described by the function Xj(a\b) (Eq. (27)), which is nothing but the Everett function 
associated with the Preisach description. The function Xj{a\b) represents the average value of x at which the 
generating function crosses the level h = afor the first time, starting atx = from the initial level h = b (see 
Fig. 6b). The description of hysteresis is reduced to the solution of this particular level-crossing problem for 
the stochastic process. 

Remarkably, this solution can be worked out in closed analytical form. To this end, let us start from the 
description of the process in terms of its Ito stochastic differential equation [26] 



where dWx represents the infinitesimal increment of the Wiener process W(x), A and B are independent of x 
because of the assumed homogeneity of the process, and x plays the role of time. The statistics of the process 
are fully described by the transition density P{h,x\ho,xo), giving the probability density that a sample function 
of the process take the value h at the position x, conditioned to the fact that it takes the value /zo at x = xq. The 
Fokker-Planck equation for the transition density associated with Eq. (30), P{h,x\ho) = P{h,x\ho,0), is 



dh = A{h) dx + Bih) dW^ 



X 



(30) 




(31) 
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As discussed before, the situation of interest is the one depicted in Fig. 6. The process starts, at jc = JCo = 0, 
from h = ho. We wish to determine the statistics of the x value at which the process reaches the level h = a 
for the first time, in the hmit ho b. This is obtained by solving Eq. (31) under the initial condition P(h,0\ho) 
= 8(h - ho), together with the assumption of absorbing boundary conditions ath = a and h = b, and then by 
taking the limit ho b. The mathematical details of the analysis are discussed in Appendix C. The solution 
for Xi{a\b), expressed in terms of the function 



V/^(m)= exp 



-2} 



"f A{u) 



du 



(32) 



reads 



Kr- ^ ^ 



{b,a] b 



du 

B^(u)\l/(u) 



(33) 



where 



(34) 



III.3. Preisach distribution associated with a given homogeneous process 

The quantity Xj{a\b) given by Eq. (33) coincides with the Everett function of the Preisach model 
associated with the homogeneous stochastic process. Therefore, as discussed in [9], the Preisach distribution 
ix{a,b) is given by 



(35) 
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By deriving Eq. (33), one finds 



[b,a] b 



du 



JB\u)\ir(u) 



(36) 



that is, taking into account Eq. (33), 



(37) 



[b,a] 



Let us calculate the Preisach distribution associated with some typical stochastic processes. 

A. Wiener process. The Wiener process is described by A(h) = 0, B(h) = 1 (see Eq. (30)). Therefore, we 
obtain from Eqs. (32) and (34), 



V(m) = 1 



(38) 



By inserting these expressions into Eqs. (33), (36), or (37) we find 



ll(a,b) = 



1 



1 



(39) 



X^(a\b) = -(a-b) 



The Preisach distribution is simply a constant and all hysteresis branches are parabolic (Fig. 8). 



B. Wiener process with drift. By this, we mean the case where by Aih) = 1/2^, B{h) = 1, with ^ > 0. We have 
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V(M) = exp(-M/(^) 



(40) 



The Preisach distribution and Xj(a\b) are given by 



ll{a,b) = 



sh^ X 



X = 



a — b 



(41) 



The Preisach distribution depends on the difference {a-b) only, and tends to the value 1/3 when {a-b) 0, in 
agreement with Eq. (38). Typical hysteresis branches calculated from Eq. (41) are shown in Fig. 9. 



C. Ornstein-Uhlenbeck process. In this case, A{h) = -h/^, B{h) = 1, with ^ > 0. We find 



-b^ 



(42) 



where ^{a,c;x) is the confluent hypergeometric function. The Preisach distribution and Xj(a\b) are obtained 
by inserting these expressions into Eqs. (33), (36), or (37). 



IV. CONCLUSIONS 

The formulation developed in the previous sections is general enough to offer various possibihties for 
further studies and apphcations. From the mathematical viewpoint, the basic issue is the role of return-point 
memory in the functional integration approach developed in Section 11. We know that return-point memory 
is inherent in the formulation, and it is natural to ask under what additional conditions, if any, an arbitrary 
scalar hysteretic system exhibiting return-point memory can be described through Eqs. (10)-(13), by 
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choosing appropriately the space G and the measure n. For the moment, we do not have a general answer to 
this basic question. 

From the physical viewpoint, the results obtained in the case where the measure is associated with a 
stochastic process are of direct interest to all those situations where some dominant degree of freedom, say X, 
evolves in a random free energy landscape, and the associated dynamics are overdamped. By this we mean 
that X obeys an equation of the form 



dX dF{x) 

r-- = H{t) ^ (43) 

dt dX 



where H(t) is the time-dependent driving field, F(X) is the free energy of the system and y> is some typical 
friction constant. Under small enough field rates, the solutions of Eq. (43) - once expressed in terms of //as a 
function of X - precisely approach the behavior shown in Fig. 1 [12,29], so that our formulation can be 
directly apphed, if one knows the statistical properties of the free energy gradient d F/dX. A particularly 
important example is the motion of magnetic domain walls in ferromagnets, where Eq. (43) often provides a 
good physical description, and various forms of structural disorder (point defects, dislocations, grain 
boundaries, etc.) are responsible for the random character of d F/dX. There are a series of classical papers in 
the literature [30,31], where the domain wall picture has been apphed to the prediction of coercivity and 
magnetization curve shapes, starting from some assumption about the properties of F{X). Equations (33) and 
(36) provide a general solution for the case where the process dF/dX is Markovian, continuous, and 
homogeneous, hi particular, the proven equivalence of Markovian disorder to the Preisach model gives a 
sound statistical interpretation of the latter in terms of stochastic dynamics in quenched-in disorder. In this 
respect, the extension of the analysis of Section m to non-Markovian and/or non-homogeneous processes 
would be of definite interest, as a way to provide quantitative predictions of hysteresis features under more 
realistic conditions and to indicate in what direction one should generalize the Preisach model in order to 
improve the macroscopic description of hysteresis generated by various forms of structural disorder. 
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APPENDIX A. PROPERTIES OF higix;xo) AND x^ih-^o) 

The function /i[g](x;xo) defined by Eq. (3) is continuous and non-decreasing with respect to x and it is 
continuous with respect to xq. Its main properties derive fi-om the following theorem. 

Theorem A.l. Given any h g [hL,hu\^ there exists at least one equilibrium state x e [xl,X[/], such that /Z[g](x;xo) 
= h. 

Proof. The theorem holds by definition for h = g{xo), h = hi, and h = hu. Then, let us consider the case g{xo) < 
h < hu- Let us introduce the set 

C/[g](/i,Xo) = {xe [xo,Xj/] : g{x) = h] (A.l) 

Because of the continuity of g, U[g\ is closed and will therefore contain its minimum x„. This is the smallest x 
e [xo,X[/] for which g{x) = h. Because g(xo) < h, then g{x) < h for any x e [xo^x^), i.e., max{g(M) ; u e [xcx^] } 
= g(x„) = h. Therefore, the equihbrium state x„ is such that /i[^](x„;xo) = h. The proof of the case hi < h < 
g(xo) is analogous. The only difference is that one must consider the maximum Xm of the set 



The function X[g](/z;xo) is the inverse of h[g]{x;xo). In fact, according to Theorem A.l, h[g]{x[g]{h;xo);xo) = h (see 
Fig. 2). It monotonically increases with h and, as a rule, it is not continuous in h. However, it is continuous in 
xo, because both max L^g] and min U[g] are continuous in xo, as a consequence of the continuity of g. Notice 
that the graph of X[g] consists uniquely of equilibrium states, that is, of points of the generating function g{x). 



Lig](h,Xo)^{x&[Xi^,XQ] : gix) = h} 



(A.2) 



Following the results of Theorem A.l, let us introduce the function X[g](h;xo) defined as follows: 




if h,<h<g(x,) 



(A.3) 
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APPENDIX B. ORDERING PROPERTIES OF T[g](/i) 

As discussed in Section 11.1, given the equilibrium states Si = (huXi) and 52 = (hi^i), we say that Si < S2 if 
xi < X2 in the usual sense. The set of equilibrium states is totally ordered with respect to this relation. Before 
considering the theorems deriving from the existence of this ordering relation, let us prove the following four 
lemmas, involving the sets U[g] and L[g] defined by Eqs.(A.l) and (A.2). 

Lemma B.l. Given xa ^ xb and hL^h< min{g(xA) , gixs)}, then max L[g](h,XA) ^ max Lig](h,XB). 

In fact, Lig]{h,XB) = L[g]{h,XA) u Lab, where Lab = {x g {xa,Xb] : g{x) = h}. Lab is empty or contains elements 

that are all greater than any element of Li^]{h,XA). In both cases, the lemma is proven. 

Lemma B.2. Given xa ^ Xb and hu>h> max{g(xA) , gixs)], then min JJ[g]{h,XA) ^ min U[g](h,XB). 

In fact, U[g](h,XA) = U[g](h,XB) u Uab, where Uab = [x g [xa,Xb) ; g(x) = h}. Uab is empty or contains elements 

that are all smaller than any element of U[g]{h,XB). hi both cases, the lemma is proven. 

Lemma B.3. Given Xa < Xb and gfe) >h> g(xB), then max L[g](h,XA) < min U[g^{h,x^. 

In fact, given any xe L[g](/z,XA) and x' e U[g](h,XB), x <Xa^Xb^ x'. This will hold in particular for x = Max 
L[g](/i,XA) and x'= min Uig]{h,XB), which proves the lemma. 

Lemma B.4. Given xa ^ Xb and g(xA) <h< g(xs), then min U[g](h,XA) ^ max L^^{h,XB). 

In fact, the set Zab - {x& [xa,xb\ : g(x) = /z} is not empty, which means that min U^g]{h,XA) - min Zab and max 
L[g]{h,XB) = max Zab- Because min Zab ^ max Zab, the lemma is proven. 

Theorem B.l. Given the equilibrium states Sa and Sb, Sa ^ Sb, then T[g](h) Sa ^ T^[g](h) Sb for any h e [hL,hu]. 
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Proof. Let us express the states as Sa = Qia = gfe), Xa) and Sb = {hs = gixs), Xb), with Xa ^ Xb- According to 
Eq. (4), the application of T[g](h) changes Xa and Xb into X[g](h;xA) and X[g](h;xB). By applying Lemmas B.l- 
B.4 to the definition of X[g]{h;xo) (Eq. (A.3)), one finds that X[g]{h;xA) ^ X[g]{h;xB) for any h e [huhu]. 

Theorem B.2. Given the equilibrium states Sa and Sb, Sa ^ Sb, and the inputs Ha and hs. Ha ^ Hb, then T[^](/za) 

Proof. Let us express the states as Sa = (Ha = gfe), Xa) and 5^ = (^b = gfe). ^s), with xa ^ x^. From Theorem 
B.l, we have that X[g](hA,XA) ^ X[^](/za;xb). On the other hand, X[g](hA,XB) ^ because of the 

monotonicity of X[g](h;xo) with respect to h. Therefore, X[g](hA;xA) ^ X[g]{hB',XB). 

Theorem B.3. Let us consider the initial equihbrium states Sa{0) and 5^(0), Sa{0) ^ Sb{0), and let us apply to 
them the two ordered input histories /za(0 ^ /zb(0- Then, at any subsequent time / > 0, 5a(0 ^ *b(0- 

Proof. Let us consider the evolution of the two states SA(t) and 5^(0 first from ? = up to the time of the first 
reversal of hA(t) or hB(t), then from this time to the time of the second reversal of /ia(0 or and so on. By 
considering that initially Sa(0) < Sb(0) and by applying Theorem B.2, we find that order is preserved in the 
first interval and that the states at the end of the interval are still ordered. This permits one to conclude that 
order is preserved also in the second interval, and so on. 



APPENDIX C. SOLUTION OF LEVEL-CROSSING PROBLEM 

The probabihty current associated with Eq. (31) is: 



j{h,x\hQ)=A(h)P{h,x\hJ — — [b ' (/j)P(/j, x I /jq )] (C. 1) 

2 dh 
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The rate at which the process leaves the interval [b,a] starting from h = hoatx = Ois obtained by integrating 
Eq. (31) over h. One obtains 




(C.2) 



Equation (C.2) shows that the probability current at the boundaries is just proportional to the probability 
density that a level-crossing event takes place at the position x. As mentioned before, we are interested in 
level-crossing through the upper boundary h = a, described by the probability current J(a,x\ho). According to 
Eq. (C.l), the functional dependence of J(a,x\ho) on x and ho is the same as that of P(a,x\ho) = P(a,x\ho,0) = 
P{a,0\ho,-x). This means that J{a,x\ho) obeys the backward equation [24] 



The total probabihty, %,a](/?o), that the process leave the interval [b,a] through h = a is given by the 
expression 



By integrating Eq.(C.3) from x = to x = 0°, and by taking into account that J(a,0\ho) = J(a,oo\ho) = 0, one 
finds that %,a](/io) satisfies the differential equation 



d d 1 

— y(a,;c I /Iq)- A(/jJ— 7(a,;c I /jJ--5^(/iJ— yy(a,;c I /io)= 



(C.3) 




(C.4) 








+ A{h,) 



dh. 



^0 



= 



(C.5) 



with the boundary conditions 
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(C.6) 



deriving from the fact that the process will certainly cross the boundary h-aifit starts from that same level, 
that is, from ho — > a, whereas it will never reach h = a if it starts from ho — > b, because in that case it will 
certainly cross the boundary h-b first. The solution of Eq. (C.5) is then 

^ib,a] (K ) = - du (C.7) 

'^{b,a] b 

where i/A^w) and Kih,a\ are given by Eq. (32) and Eq.(34), respectively. 

The probabiUty density p[b,a]{x\ho) that the process reach the boundary h = asA. the position x is given by 



/ , , X j{a,x \ hr.) 

^[bAK) 



and the mean value X[b,a](ho) of the level-crossing position is 

^[b,a] (^0 ) = j " Pib,a] {u\K) du = \7-\ \u j[a,u\hQ) du (C.9) 

^\b,a] V"0 / 

By definition, the conditional probability density T{a,x I b,Q) of Eq. (27) is given by the limit of Eq.(C.8) for 
ho — > b, that is 

T{a,x\b,Q))=p^,^^{x\b) (C.IO) 
and the function Xj{a\b) defined by Eq. (27) is accordingly given by 
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Xj.(a\b) = x^i,^^(b) 



(C.ll) 



By multiplying Eq.(C.3) by x, by integrating it from x = to x = 0°, by taking into account that J(a,0\ho) 
J{a,oo\ho) = 0, and by making use of Eq. (C.4), one finds that the function 



(C.12) 



obeys the differential equation 



dh^ dhg 



(C.13) 



with the boundary conditions 



(C.14) 



deriving from the fact that X[b,a](a) = by definition, whereas %,a](i>) = because of Eq.(C.6). Equation 
(C.13) is a linear, non-homogeneous first-order differential equation for df[h,a]/dho. The solution reads 



w{K\c-2\^^^du 



B\u)\if{u) 



(C.15) 



where and \j/{u) are given by Eq.(C.7) and Eq. (32), respectively. The constant of integration C can 

be expressed as 
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c = 



'df, 



[h,a] 



dhn 



(C.16) 



By deriving Eq.(C.12) with respect to ho, by taking into account that 7C[i,^a](b) = and [d7i[ba]/dho\b 
yA,b)/K[i,,a], and by making use of Eq.(C. 1 1), we obtain 



lb,a] 



y/(b) y/(b) 

^ib,a](b) = — X^(a\b) 



{b,a] 



(C.17) 



Equations (C.16) and (C.17) permit one to write Eq.(C.15) in the form 



dhn 



1 ^[.,«] iB\u)xif{u) 



du 



(C.18) 



By integrating Eq.(C.18) from b to ho and by inverting the order of integration in the double integral, one 
obtains 



(^0 ) = (« I ^K,a] (^0 )- j \y{u)du [ \lf{u)du 



\.b,a] b 



du 



(C.19) 



where use has been made of Eq.(C.7). Taking into account that /[i,_a](a) = and 7C[b,a](a) = 1, one concludes 
that Xj(a\b) must be equal to Eq. (33). 
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Fig. 1 - Free energy F(X) with multiple minima and corresponding gradient dF/dX. The dashed line 
represents the hysteretic behavior of H(X) obtained from the stability condition H = dF/dK. 




Fig. 2 - Generating function g{x) with illustration of the envelope construction (broken line) associated with 
the function h[g\{x;xo) and its inverse X[g](/z;xo). 
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Fig. 3 - Action of evolution operators T[^](/i) associated with a given sequence of input reversals. Initial state 
is lower saturation sl - {hi^i). Reversals take place at s\ - T[g](/zi) sl - {hi^ig^ihuXi)), S2 - T[g](/z2) s\ - 
(/z2,X[g](/z2;xi)), and so on. 




Fig. 4 - Illustration of intervals with different kinds of stability. AB: weakly stable; A'B': strongly stable; 
CD: weakly unstable; CD': strongly unstable. 
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Fig. 5 - Typical generating function associated with the Preisach model (solid line), and corresponding 
envelope construction (dashed line, see also Fig. 3). The dotted lines give examples of different weakly 
unstable behaviors giving rise to the same hysteresis properties. 




Fig. 6 - Left: example of stochastic process sample functions involved in the study of level-crossing through 
the boundary h - a (thick line) ox h- b (thin line), starting from the initial condition x - xq &t h - ho. Right: 
same as before, in the particular case where b and crossing through h- a only is considered. 
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Fig. 7 - Level-crossing problems to be solved to calculate hysteresis in Markovian process. Top (field 
increasing from hi toward h{): level-crossing must be considered in the interval [hi,h,], with known 
distribution Pl(xl) at the lower h boundary hi. Bottom (field decreases from hi toward h2, not shown): level- 
crossing must be considered in the interval [h„hi], with known distribution pi{xuh{) at the upper h boundary 
hr. 
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Fig. 8 - Typical hysteresis curves for Wiener process, calculated from Eq. (39). 




Fig. 9 - Typical hysteresis curves for Wiener process with drift, calculated from Eq. (41). 
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